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Abstract 

In addition to the well established large atmospheric angle, a large solar angle is probably present 
in the leptonic sector. In the context of the see-saw and by means of a bottom-up approach, we explore 
which patterns for the Dirac and Majorana right-handed mass matrices provide two large mixings in 
a robust way and with the minimal amount of tuning. Three favourite patterns emerge, which have a 
suggestive physical interpretation in terms of the role played by right-handed neutrinos: in both solar 
and atmospheric sectors, either a single or a pseudo-Dirac pair of right-handed neutrinos dominates. 
Each pattern gives rise to specific relations among the neutrino mixing angles and mass differences, 
which lead to testable constraints on U e %. The connection with the rate of LFV charged lepton decays 
is also addressed. 
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1 Introduction 



Combined results ||] from SNO and Super-Kamiokande || suggest that another large mixing angle 
is present in the leptonic sector, in addition to the well established large mixing angle involved in 
atmospheric oscillations 0. CHOOZ data J5|, on the contrary, constrain the third angle #13 to be 
smaller than the Cabibbo angle. Such an interestingly different pattern of lepton mixings as compared 
to quark mixings should be interpreted as the result of some specific property characterizing the physics 
beyond the Standard Model. 

The see-saw mechanism Q has emerged as the most elegant explanation of the smallness of neutrino 
masses: Ai v = Y T M^Yv 2 , Yv and M being respectively the Dirac and the Majorana mass matrices 0. 
Triggered by such impressive experimental results, one would like to understand if the requirement of a 
large solar angle leads to some specific realization of the see-saw mechanisms. Then, one could look for 
theoretical implications - like those on the underlying flavour symmetry - and also for phenomenological 
implications - like the magnitude of U e 3 and the rate of lepton flavour violating (LFV) charged lepton 
decays g, g. Many see-saw neutrino mass models have been proposed in the last years (for a review 
and a set of references see e.g. H), in particular in connection with flavour symmetries and grand 
unified theories. From the model building point of view, with an hierarchy between m@ = \J Am 2 tm 

and m.Q = yj Am|, it is much more challenging to end up with a large mixing than with a small one. 
This holds for the atmospheric angle, so that it holds a fortiori when requiring a large solar angle 



too. Indeed, either one imposes various tunings between the many elements of Y and M [10], or one 
has to deal with apparently intricate patterns for Y and M which could only be obtained from highly 
non-trivial flavour symmetries p[ |i~2| , 13, 14, 15]. 



To overcome the unavoidably fragmentary picture which arises by considering the many specific 
models 0, one should investigate in a model independent way which essential features of Y and M 
can account for two large angles and two mass scales. In this paper we tackle this problem from a 
bottom-up point of view and we identify those Y and M patterns which give two large mixings because 
of their own structure with the minimal amount of tuning [J Since their structure is just the reflection 
of the different roles played by right-handed neutrinos, from the analysis of such 'natural' patterns 
it will turn out that few suggestive possibilities for such roles. As it will be discussed, interesting 
phenomenological implications also follow, in particular for the magnitude of U e %. 

At low energy and in the basis were charged leptons are diagonal, the neutrino sector is described 
in terms of nine physical quantities, while at high energy, that is above the scale where right handed 
neutrinos are integrated out, it is specified by eighteen. Casas and Ibarra [18] identify the hidden 



parameters of the see-saw with the six parameters of a complex orthogonal matrix acting on the left 
side of L = M~ l l 2 Yv and the three right handed neutrino mass eigenvalues. Of course, it is unrealistic 



at present to follow a complete bottom-up approach to reconstruct Y and M [19]. In fact, even if 



other observable quantities, like the rate of LFV processes Hlq , 2C, 21], are sensitive to the hidden 
see-saw parameters (see e. g. [|jj|, |22|), such observables - even if measured - provide upper limits, so 
that one is still far from having any direct access to the elements of Y and M. 

In order to fill this gap between the number of observables at low and high energy, one needs 



1 We are dealing here with the simplest version of the see-saw mechanism, obtained by just adding three right-handed 
neutrinos to the matter content of the Standard Model. 
2 For previous interesting attempts see J^J. 

3 We will not consider the extended flavour democracy scenario (see e.g. [[nj and references therein). 
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some 'theoretical salt'. The construction of specific models [^] according to the rules dictated by 
some flavour and/or GUT symmetry is a top-down approach to the problem of neutrino masses and 
mixings and the only ambiguities are related to the predictive power of the theory. In this sense, 
the neutrino oscillation data can test the different theoretical high energy models, but of course they 
cannot select a particular one. 

In the bottom-up approach of this paper, we adopt instead the requirements of robustness and 
economy as additional ingredients. Indeed, our aim is to isolate those {Y, M} which account for 
neutrino mixings and the hierarchy m@ >> ijIq in a natural way. In this spirit, we look for the 
structures of Y and M which are stable, in the sense that small perturbations in the parameters 
just induce small corrections in the physical observables. Thus, a large mixing angle cannot result 
from a tuning of the many elements of Y and M, which would be presumably unstable, but rather 
from a simple relation between a couple of matrix elements. As another consequence, we exclude the 
nearly degenerate neutrino mass spectrum because the mixing angles are very unstable under radiative 
corrections p4| ]. The approach is also economical as it minimizes the number of the conditions on 
the parameters, while allowing most of the parameters to remain generic. One can derive some 
generic predictions in the framework of each pattern in the form of quantitative constraints between 
experimental observables. In this sense, we may decide if one generic pattern is plausible or disfavoured 
by present and future data. We carry out our analysis for the two viable scenarios for the neutrino 
mass spectrum, namely the hierarchical (Hi) and inverse hierarchical (iH) case. 

We first apply the stability criteria to the analysis of the atmospheric angle. In Section 2, we 
display the patterns of M v consistent with large atmospheric mixing and m@ >> tuq for these two 
scenarios. Then, we discuss how to obtain them with the least amount of tuning. In the (Hi) case, the 
naturalness criteria for the atmospheric sector are satisfied by the well known dominance hypothesis 



25. 12, 13]: the mass scale of the heaviest light neutrino is essentially provided by only one right handed 
neutrino, and the fact that the angle is large is obtained by imposing that its Yukawa couplings to 
the /i and r have comparable strength. The dominance hypothesis has interesting phenomenological 



21, p| and theoretical [13] consequences. Also for the (iH) case, as will be shown, it is possible to 



recognize a dominance mechanism which would be solidly realized by a pseudo-Dirac pair of right- 
handed neutrinos having a suitable hierarchy in their Yukawa couplings. Models for (iH) have been 



studied in particular in association with the breaking of L e — — L T [£7], 28 



Along this line of reasoning, the next step is to investigate which pattern of the sub-leading 
elements of Ai u accounts for the presence of a large or maximal solar angle in the most natural 
way. This is done in Section 3. Dealing with sub-leading contributions is a delicate task in the 
presence of several small parameters. Moreover, one has to take into account the CP violating phases 
|3l| ] potentially present in Y. We carry out this program in a quantitative and systematic way. Of 
course some results are already known as empirical rules from many previous analyses, in particular 



from the model building point of view (see, e.g. |25|, 11, 12, [16|]). We perform the calculation in a 
systematic way to see explicitly how much these intuitive ideas are reliable and if there are important 
exceptions. Fortunately, a more quantitative reappraisal is already possible with reliable approximate 
expressions, where to recognize immediately the leading and the sub-leading contributions, keeping 
under control the order of magnitude of the negligible sub-subleading contributions. These expressions 
show that a large solar angle represents a strong constraint and suggest peculiar structures for Y and 
M. This also explains why it is non trivial to build models which account for two large angles. It 
turns out that the quantitative analysis supports patterns which have a remarkably simple physical 
interpretation in terms of how right-handed neutrinos couple to each light neutrino flavour eigenstate. 
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We then describe in some detail these suggestive patterns. In the (Hi) case, as we have already 
mentioned, the atmospheric mass scale and mixing angle essentially arise from just one of the right- 
handed companions, so that the solar sector is directly linked to the properties of the other two 
right-handed neutrinos. Two basic possibilities automatically lead to a sufficiently large solar angle 
which we will refer to as double dominance and pseudo-Dirac for short. The first case presents also in 
the solar sector a mechanism analogous to the one at work for the atmospheric sector: one among the 
remaining two right-handed neutrinos is essentially responsible for the solar mass scale and the large 
solar angle is obtained when it couples with the same strength to the e and to the \x and r combination 
orthogonal to the heaviest neutrino eigenstate. To approach a maximal mixing would need a precise 
equality of the couplings. In the double dominance each light neutrino has its preferred right-handed 
companion. The second possibility for having a stable large solar angle, arises when the solar sector 
of the light neutrino mass matrix possesses a pseudo-Dirac structure. The solar angle approaches its 
maximal value when the Dirac structure is stressed. This naturally arises when the two remaining 
right-handed neutrinos possess a pseudo-Dirac structure and a suitable hierarchy of Yukawa couplings. 
In the (iH) case, the mass of the pseudo-Dirac pair is close to the atmospheric scale; m & being much 
smaller, the pseudo-Dirac structure has to be well approached. As a consequence, the deviation of 
the solar angle from tt/4 has to be very tiny. The relation is so strong that only the LOW solution is 
really viable. 

From the point of view of introducing the minimal amount of tuning, these three scenarios are 
preferred. They are all characterized by the dominance of either a single or a pseudo-Dirac pair of 
right-handed neutrinos in both the solar and atmospheric sectors. On the contrary, the simplest flavour 
symmetry, a plain U(l) symmetry [^] with all charge assignments of the same sign, would predict a 
democratic dominance pattern, with all right handed neutrinos contributing roughly the same to all 
light neutrino masses. Then, the above considerations show that the most natural structures from a 
bottom-up perspective are realized only by means of non-trivial flavour symmetries such as t/(l)'s with 
holomorphic zeros from superymmetry O, 14] or non-abelian flavour symmetries [15]. Alternatively, 



one could consider the above three possibilities as recipes for model building. Of course, deviations 
from these structure require introducing more correlations among the parameters. 

The above natural patterns for Y and M are also interesting from a phenomenological point of 
view. Being quite constrained, it is straightforward to look for their predictions. In particular, those 
for C/e3 are remarkably different: in the (Hi) case, \J e % could not escape detection if the experimental 
sensitivy reaches few percent, while in the (iH) case, on the contrary, only LOW is viable and U e z is 
predicted to be smaller than one per mille. All the three scenarios considered predict a rate for f3(3Qi> 
too small to be observed. In Section 4 we discuss the connections with the LFV decays /i — ► e^f and 
r — ► /i7, emphasising the strong dependence on the absolute scales of right-handed neutrino masses 
and on the other hidden see-saw parameters. On the contrary, sizeable electron and muon electric 
dipoles would provide an indication for a hierarchy in the right-handed neutrino spectrum [p3|. 



In section 5 we show that the hidden parameters identified with the complex orthogonal matrix 
admit an interesting physical interpretation as a dominance matrix. Outlook and conclusions are 
presented in Section 6. 
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2 Accounting for the Atmospheric Neutrino Mixing Oscillations 



2.1 The See-Saw Parameters 



In the see-saw model || the effective neutrino mass matrix is given in terms of the Dirac mass matrix 
Yv and the Majorana mass matrix M for right-handed neutrinos as 

mi \ 2 

m 3 / 

where y is defined in the basis where the charged leptons are the mass eigenstates (e, /i, r) and, for 
definiteness, we choose the basis of M eigenstates (M x , M y , M z > 0) for the right-handed neutrinos. 
Their ordering will be defined below. The complete determination of Ai u would provide us with 
9 parameters: 3 eigenvalues of M v , 3 angles and 3 phases of U. Instead, the see-saw model has 
21 parameters: 18 complex numbers of the matrix Y and 3 eigenvalues of M. At the lepton flavour 
violation scale, in the framework of the see-saw model (i.e., the SM or the MSSM plus the right-handed 
neutrino masses), 3 phases in Y can be reabsorbed in the definition of the leptonic fields. There are 
18-9=9 more parameters in the see-saw model than those to be measured in the matrix ]\A. V . 

Let us define the matrix 

L = M~ 1/2 Yv = 

such that M u = L T L and the bases are defined as above (3 phases in the 9 complex matrix elements 
are unphysical at low energies). The see-saw mechanism involves L and the M eigenvalues. If CP 
is conserved in the lepton sector of the whole theory (i.e. at energies above the scale where right- 
handed neutrinos decouple), the phases of the fields can be chosen to have Y, M and A4 U to be real 
matrices. Negative M (M u ) eigenvalues correspond to CP odd right-handed (left-handed) neutrino 
eigenstates. It is worth stressing the following property: M and A4 U possess the same number of 
negative eigenvalues. Actually, we shall prefer here the convention where M and M v eigenvalues are 
all positive, and multiply the lines in Y and L that correspond to a CP odd right-handed neutrino 
state by a phase i. In general however, the phases may be important and the L matrix elements are 
complex, while M x , M y , M z are positive. 




As noticed by Casas and Ibarra |18[| , extracting L from Jv[ v suffers from the ambiguity 

M u = L T L = {RL) T {RL) , (3) 
where R is any complex orthogonal matrix, i.e., any element of 0(3, C). More explicitly, one has 

Mvij = XiXj + yiyj + ZiZj, (4) 



which is clearly invariant under 







f Xi 


y% 




y% 


Zi ) 




V Zi 



1,2,3, (5) 
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where R can be parameterized as 

R = e aJ 'e^ J "e lfJ ', (6) 

with a, (3, 7 three complex numbers and J x ,y,z are the 0(3) generators. Therefore the solutions L of 
eq. (H|) can be written as 

/ x 3 \ 

L = R y 2 y 3 , (7) 

\ zi z 2 h ) 

where R £ 0(3, C) remains arbitrary and the matrix elements x 3 , ...,-?3 (among which 3 are real, the 
other 3 being complex) can be determined from M u , up to sign ambiguities. Of course, both R and 
the eigenvalues of M are physical in the framework of the see-saw model, but other experiments are 
needed to access them, as we also discuss below. 

In some cases discussed below, the matrix R defined by Casas and Ibarra (which should not be 
confused with R) can be interpreted as a dominance matrix, which associates each light neutrino 
mass eigenvalue nii with a given combination of right-handed neutrino masses. In general however, 
its interpretation is not transparent because it is not unitary. For this reason, we prefer to work with 
(|) rather than with (0) and ©. 



2.2 See-Saw Realizations of Mass Hierarchies and Large Mixings 

Our aim in this section is to identify the patterns for the see-saw parameters - more exactly, the 
patterns for L, since we do not assume any a priori hierarchy for the right-handed neutrino masses - 
that fit in the most natural way the experimental results. Because of the surplus in the parameters, 
one must adopt a simplifying assumption and it looks reasonable to assume each feature of Ai u to 
be due to a simple relationship in L, such as two parameters being almost equal or very different, 
rather than to assume an elaborate conspiracy between the many elements of a matrix. We also apply 
the same considerations in linking the properties of L to those of Y and M. We first sketch the two 
patterns that allow to minimize the tuning of the seesaw parameters. 

The first feature is the hierarchy between m% and nvL. A way to account for one eigenvalue being 



much larger than another, is to assume a dominance mechanism [25, 12, 13 1, in which one right-handed 
neutrino, associated to a line in L which we take to be z, gives the dominant contribution to the larger 
eigenvalue m,j of M. v . Then 



J2((LU)jif™ «E(^) 2 ^ 2 , (8) 



and (phases are unimportant here) the corresponding eigenstate is approximately ZjU 3 /z, so that 

U e i^zi/z, U fli ~z 2 /z, U Ti ~z 3 /z . (9) 

Within the dominance assumption, M u ij is given by ZiZj plus corrections of the order of the smaller 
eigenvalues. Whenever some Zi/z becomes small as compared to these corrections the corresponding 
expression in (|9|) has to be corrected. Thus, the larger eigenvalue of M. v is here due to a right-handed 
neutrino with some larger Yukawa couplings or a relatively small mass. 
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The second feature of neutrino oscillations is the favoured large mixing angle solar solutions. A 
large mixing angle is naturally obtained with a pseudo-Dirac structure of Ai u , i.e. when the dominant 
entries are off-diagonal (as discussed below, this cannot be applied to the atmospheric mixing angle if 
almost degenerate neutrinos are excluded). In this case at least two mass eigenstates are very close in 
mass, with CP-violating Majorana phases differing by tt/2 (opposite CP parities in the CP-conserving 
case), and strongly mixed. Such a structure requires strong correlations between the couplings of at 
least two right-handed neutrinos. The most natural origin for these correlations is the existence of a 
pseudo-Dirac pair of right-handed neutrinos - recall that, in the CP-conserving case, there should be 
as many right-handed neutrinos with negative CP parities as light neutrinos with negative CP parities. 
We concentrate on the 2x2 pseudo-Dirac sector of the mass matrices. Then there is a basis in which 
M takes the form: 

jQ- (+ sman corrections) , (10) 

where we have chosen the pseudo-Dirac pair of right-handed neutrinos to correspond to two lines in 
L, e.g., x and y. Neglecting the smaller contributions (deviations from the pseudo-Dirac pattern) one 
gets for the L matrix elements: 

{Y 2i + Y u ) i(Y 2i -Y u ) 
Vi = 1= — v , Xi = 7 == — v , (11) 

where the Yukawa matrix Y is in the basis defined by (|l0| ) . Eq. (|Tl| ) leads to a pseudo-Dirac structure 
for Ai u , with mi ~ — m 2 , provided the hierarchy among Y couplings is as follows: 

Y 21 » Y u and Y 12 » Y 22 (or Y 2j ^ Y Xj ) . (12) 

so that \y\ + ix\\ <C \y 2 + ix 2 \ and \y 2 — ix 2 \ <C \y\ — ix\\ or vice versa. 

The diagonalization of the resulting M u yields a large mixing angle which becomes closer to 
7r/4 as the pseudo-Dirac configuration is better approached. Therefore, a close to maximal mixing 
angle in neutrino oscillations (as indicated by the solar data in the LOW region) suggests a pseudo- 
Dirac structure in M and in M u , which requires in turn that some particular Yukawa couplings are 
predominant. However this option cannot be advocated for the atmospheric oscillations, because it is 
not possible to accommodate solar neutrino oscillations within a 3- neutrino scheme with m 2 ~ — m.3 > 
m@. 

As already stressed, an almost degenerate mass spectrum - perhaps the most obvious explanation 
for the large mixing angles - is not yet excluded by experiments. But it is generically disfavoured be- 
cause the large mixing angles become strongly scale dependent, unless some compensation mechanism 
is introduced 



2.3 Patterns for M v 

Once solutions with three almost degenerate eigenvalues are excluded as unstable under the RGE 
evolution, the largest mi eigenvalue must be 0{m@). Furthermore, the matrix M. u must be consistent 
with the experimental facts: 

a) mQ = m\ — m\ <C mj|, namely the hierarchy between solar and atmospheric neutrino mass 
differences; 

b) large atmospheric mixing angle #23 ~ vr/4 and good evidence for large solar mixing angle 0\ 2 ~ 7r/4 
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as well; 

c) the CHOOZ constraint |J7 e3 | < 0.2. 

There are basically two patterns for A4 U that are consistent with these properties (the choice of 
the ordering i = 1,2, 3 of the nii is suitable for the standard conventions), which we now discuss in 
turn. 



2.3.1 Hierarchical pattern: 7713 ~ ?ti@ 3> mi, mi 

The properties a), b), c) lead to a pattern for M u consistent with the assumption of the dominance 
of the atmospheric oscillations by one right-handed neutrino, namely (the intervals at 90% C.L. are 
taken from Ref. §4|), 

M vij = ZiZj(l + 0(p)) (13) 

with 

3 

^2 z f - m ® = ( 4 - 8 ) x l(r 2 eV 
i=i 

— ~ tan 023 = (0.6- 1.7) (14) 
\z 1 /z 3 \ < tan < 0.2 

and 0{p) denotes a matrix whose entries are at most of order p, with p = 7712/7713. Thus this pattern 
is characterized by dominant matrix elements of order 77i@ in the (2,3) submatrix of M u , so that 
Z2 ~ z 3 ~ \J m@/2, while the other are smaller since z\ < \fm©,/2tan6i 3 . This means that one 
right-handed neutrino dominates the atmospheric oscillations. According to the values of Xi,yi, there 
are several possibilities to obtain m and 612 from the other two right-handed neutrinos, as discussed 
in the next section. 

In a bottom-up approach, it is convenient to redefine L by a rotation 

(x 2 , x 3 ) — > (x-, x + ) = (cos6> X2 - sin 6 x 3 , cos 6*0X3 + sin #0^2) , (15) 

where O = arctan^/^), and analogously for (2/2,2/3) (y-,y+), (z 2 ,z 3 ) -> (^_,z + ). Since 9 ^ 6*23, 
the (+) components are approximately aligned with the heaviest neutrino mass eigenstate z/3, while 
the (— ) components are orthogonal. After this rotation one has z_ = and z+ ~ y/m@, which implies 

M v+ -, M v — <S M v++ ~ Til® . (16) 

The atmospheric angle 6*23 differs from 9q by a small angle that can be treated perturbatively as we 
do in the next section. In the rotated basis, 

' x\ 

L=\yi y- y+ I , (17) 




with z + ~ i/rn@, and all other entries are smaller than Jm^: x±,y± <C ^ l /rn@ expresses the dominance 
condition, while x\,y\,zi <^ ^/m@ reflects the smallness of tan 0i3. How small they should be depends 
on the favoured solar neutrino solution and on the value of tan #13. 

This pattern which is quite suggestive of the properties a), b), c) is not obtained in straightforward 
applications of abelian flavour symmetries to the see-saw mechanism []32]]. But it can be implemented 
through the use of holomorphic zeros or non-abelian flavour symmetries 14, 15 1. 
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2.3.2 Inverted hierarchy pattern: mi,m,2 — m@ 3> 1TI3 

Since mj — m\ = m|, <C mjg, mi ~ ±m2, but only mi ~ — m2 is consistent with the stability of 
the large mixing angle #12 under radiative corrections. We shall therefore focus on the case where 
the mass eigenstates u\ and v 2 form an approximate pseudo-Dirac pair. In this case the structure of 
the light neutrino mass matrix indicated by the data - leaving aside the now disfavoured small angle 
MSW solution - is the following: 



A-tan 2 6»i2 m 3 2 m Q \ 

Muii ~ rn@ max , U e3 , , 

V 2 m@ mi/ 



(I -tan 2 6>i2 TT m 3 m% 

M u22 33, 23 ~ m @ max 7, ) Ue3, , — 7T . , 

V 2 rn@ m| 1 

M vl3 ~ 7W^i2 ~ m@ , - ~ tan ^23 • (18) 

Ali/ 12 

Thus ^^12 and M u i3 essentially determine the atmospheric neutrino parameters. The other matrix 
elements are related to the solar neutrino parameters and are much smaller. 

Let us identify the patterns of L that are compatible with those requirements. The more natural 
way to implement the mass hierarchy between the pseudo-Dirac pair and u 3 is to associate the first 
to two right handed neutrinos (in the sense that m@ it is dominated by the z% and fji) while 7713 is 
dominated by the X{. Neglecting the contribution of the X{ to the entries of M, v , Eqs. (|l8|) can be 
rewritten as: 

Z1Z3 + VW3 - m@ sin 6*23 , 
z\z 2 + yiy2 ^ -m@cos6>23 , 

\ziZj + yiyj\ < m@ for i = j = 1, or i,j / 1 . (19) 

and also \xiXj\ <C m@ . This is precisely the pattern obtained in Section (fT^) where it is shown that 
the simplest realization of a pseudo-Dirac mass structure ([19]) is when M y and M z form a pseudo-Dirac 
pair with some hierarchy in the Yukawa couplings. 

If CP is an exact symmetry, we can choose the lepton phases such that M u is a real matrix so 
that each line in L must have matrix elements that are all real or all purely imaginary. Therefore the 
Zi (i = 1, 2, 3) can be taken to be real and the yi (i = 1, 2, 3) to be imaginary in order to generate the 
opposite signs of the masses mi and ?tt,2 , also in correspondence with the opposite CP parities of M y 
and M z (the effects of CP violating phases will be discussed in more detail later). 



In this inverse hierarchy scenario, it is useful to change to the new variables: 

Zi + iyi Zi - iyi 



(20) 



V2 V2 
so that Mvij = UiVj + UjVi + XiXj . The condition (|l2|) reads: 

Vl > 1>2,3 ^2,3 > ui ■ (21) 
As for the hierarchical spectrum one defines ± components by the rotation 

(x2, X3) — ► (x + , X-) = (cos 6 x 2 - sin 6 X3, cos 60X3 + sin ^0^2) , (22) 
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with 9q = arctan(— U3/U2). After this rotation one has ti_ = and v±u + ~ m@, which implies 

M v +-, Mui- < M u i+ ~ viu + ~ m@ . (23) 

The atmospheric angle #23 differs from #0 by a small angle that can be treated perturbatively as we 
do in the next section. In the rotated basis, 

( X\ X + X- \ 

(24) 





( Xx 


x + 


X- 






u + 














The symmetry of the solutions under yi — > — yi corresponds here to U{ <-> 



There are three possibilities consistent with (p3|): 1) \v\\ > y/m@, > \u + \ , 2) \u+\ > y / m@ > |ui| , 
3) |tt_|_ I ~ I vi| ~ y/m@. Since these three kinds of inverted hierarchy solutions correspond to the same 
pattern for A4 U , they should be connected by orthogonal matrices R. If we use the parameterization 
(0) with R = 1, one obtains from the conditions (|l9|) the solutions of type 2). A boost, i.e. a rotation 
by an imaginary angle iX in the (y, z) sector transforms the variables (u, v) as 

u-i — > e x Ui Vi —> e~ x Vi u^j —> UiVj . (25) 

When A is large enough, a solution of type 1) is transformed into a solution of type 2). In between, 
the solutions 3) can be obtained. These orthogonal transformations (]25| ) parameterize a continuum of 
equivalent mathematical solutions with quite different physics contents (see Section |j) . Notice that in 
each case M is defined to be diagonal and does not commute with R. 

1 11 

The elements of the matrix L are naturally bounded by 0(v/M i ), as far as the Yukawa couplings 
are perturbative at the scale of lepton number violation. Therefore, if e.g. some yi ~ y/m@, one has 
My v 2 /m@ ~ 5 x 10 14 GeV, and analogously for the Xi and Zj. In the hierarchical case, one has 
(with our convention that m@ is dominated by M z ) M z < 5 x 10 14 GeV, while M x and M y remain free. 
In the inverted hierarchy case, there are two possibilities: for the case 3), M y and M z are bounded by 
5 x 10 14 GeV, while in cases 1) and 2) the upper bound is lower since either v\ or u+ is larger than 
y/m@. 

Now that the large scale m@ has been identified, the oscillation parameters 6*13, #12, 623 — #0 and 
rriQ can be determined by an expansion in these small parameters that are constrained by experimental 
data in turn. This is discussed in the next Section, where CP violating phases are also taken into 
account. 



3 Origin of The Large Solar Angle 

As discussed in the previous Section, only two patterns for ]V[ V are consistent with the requirement 
of a large atmospheric mixing and m@ >> 7Bq. Along this line of reasoning, the next step is to 
investigate whether the requirement of a large - or possibly maximal - solar angle can give further 
informations on the sub-leading terms of those two patterns. It is important to take into account the 
CP violating phases potentially present in Y and so in L . In the following we carry out this program 
in a quantitative and systematic way. Of course some results are already known as empirical rules 
from many previous analyses. Fortunately, a more quantitative reappraisal is already possible with 
reliable approximate expressions. 
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In this section we adopt the following notation: 



Am B A = B (l + 0(U* 3 ,m Q /m @ ))) , (26) 

which relates the neglected contributions to the physically small parameters. 

We are seeking for sufficiently well approximated expressions for the observables contained in J^A. V ^ 
in particular for the three mixing angles and the CP violating phases of the MNS matrix U, which 
can be written in the form 

U = e ia W R(9 23 )R(9 13 ) diag(l, e"* 1)R(6 12 ) V . (27) 

where W and V are diagonal matrices of SU(3). However, in the interaction basis where charged 
lepton masses are diagonal, there is still the freedom of phase redefinition for each one of the leptons 
e, r; at the level of the matrix L this amounts to the freedom of multiplying each column by an 
arbitrary phase. 

Therefore a and the phases in W are not physical as they can be absorbed in the e, /x, r (on 
the contrary, the two Majorana phases in V and the CKM-like phase %, often denoted by 5 in the 
literature, are physical; in particular x could be measured in oscillations at a neutrino factory). Once 
this is done, it is convenient to define the matrix 

e = LR(e 23 )R(e 13 ) . (28) 

Indeed the effective light neutrino mass matrix reduces, in the basis defined by £, to 

m = £ T l = ( mi2 m 22 , (29) 
m 3 J 

and the solar neutrino parameters trig and 9\ 2 , as well as the CP-violating phases x an d V, can be 
obtained by diagonalizing a 2 x 2 matrix. 

The CKM-like phase x is determined by the condition 

\m 22 \ sin(a 2 — X) = \ m n\ sin(ai + x) , (30) 

and the solar angle is given by 

tan 2^i2 = , ; 2 mi ' , , - , , (31) 

1 77122 1 cos(a 2 - X) - l^n I cos(qi + x) 

where a 2 = &rg(m 22 ) — arg(mi2) and a\ = arg(mn) — arg(mi2)- Eq. (|30| ) has two solutions which 
differ by it, and correspondingly tan2#i2 can take both signs. However the two solutions are physically 
equivalent, as expected, since for < 6*12 < vr/2 they are related by exchanging mi and 7712- Thus, it is 
not restrictive to take ttiq = Am^ = m\ — m\ > 0, so that 1 777,22) > |t77h| corresponds to < ^12 < vr/4 
while 1 777.22 1 < 1 77711 1 corresponds to 7r/4 < ^12 < vr/2 ("dark side"). 



Eq. (|30D is graphically solved for x by the construction of the two triangles in Fig. 1. The other 



phase dependent quantities discussed in this section, in particular the denominator of (31), are also 
shown in the figure. The CP-conserving case corresponds to a flat triangle. One can easily see that, 
for fixed |mn| and 1 77722 1, the maximal (resp. minimal) cancelation in the denominator of Eq. (|3l|) is 
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Figure 1: The triangle denned by the condition on x eq- (p0[): AB = A'B = A" B = |mn|, = 
1 77^22 1 ■ The denominator of eq. (31) is A'D = I77122I cos(a2 — x) ~ \ m u\ cos(ai + x) > A"D, where 
A" D = A"'D = 1 77122! — l m n|- Thus, phases work against the large mixing. 



obtained in the CP-conserving case with 7711177122 > (resp. 7711177122 < 0). Thus the general effect of 
the CP-violating phases is to reduce the solar mixing angle with respect to the CP-conserving case 

with 7711177122 > 0. 

As for ttiq, one has one phase independent relation, 

TTlg COS 2^12 = |?7i22| 2 — |t71ii| 2 (32) 

and another one is 

77i Q sin 2#i2 = 2|?7ii2| (I77122I cos(q 2 - x) + l^nl cos(«i + x)) , (33) 

which basically gives 77i Q for near maximal mixing. 

Two effects can conspire to give a large mixing angle: |mi2| > \ m 22\ and some amount of 

cancelation in the denominator of Eq. (|3l|), with 1 777.22 1 cos(ct2 — x) ~ l m n| cos(ai + %)■ If we combine 
© and (H) into 

77i| tan 2 #i2 
1 777,12 1 2 1 1 — t an4 ^12 

one has a criterion to decide which mechanism is preferred. Indeed, if the denominator in (|3l| ) is made 
small by cancelation, the ratios in ( |34| ) are expected to be much large. Instead, if both terms in this 
denominator are small, these ratios are 0(1). This criterion will be applied to different possible seesaw 
patterns here below, by first obtaining also constraints on |r?7i2 1 - 



1 77122 1 cos(a 2 - x) + l m n| cos(qi + x) 



1 77122! cos(a 2 - X) - l m n| cos(qi + x) 



>1, 



(34) 
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3.1 Hierarchical Spectrum 



Let us first consider the case of hierarchical neutrinos. As already discussed, the associated pattern of 
My suggests that there are two dominant elements in L, which we are free to place in the third row 
of L - namely zi 3. The freedom linked to the phases W can be used to ensure that #23 and #13 are 
real. The expressions for these angles can be written as 

tan 6>23 ~ 1 9 , (35) 

z 3 z 2 

tanfl 13 « — + yw+ „ 1 + , 36 

where the variables (cc±, y±, z±) are defined in ( ^q) , y + = y + + (zi/z + )yi and y± = yi — [z\j z + )y + , 
and analogously for x+ , xi. Therefore the two phases in W are chosen to eliminate the overall phases 
in these two expressions, which are independent of the phase e lOL . 

After performing these rotations in this approximation, one gets 

X\ X- 

I = LR(9 23 )R(6 13 ) « I m y- y+ | , (37) 

x 1 x + +y l y + x-x + +y-y+ 

Z + 2+ 

where z 2 = zf ~ m@ . Then m reads 

Vi + x\ yiy- + xxx- 
m ~ I ViV— + y^L + x 2 _ I . (38) 

z 2 + y\ + x 2 + 



By replacing these matrix elements in ( 3C ) , (|3T|) , (32) and (|33|), one obtains the expressions for 
tan 2 ^i2 and thq. In the following, we shall describe separately the two mechanisms that can lead to a 



large atmospheric mixing angle: cancelation in the denominator of Eq. ( |31|) and |mi2| 3> |w-n|, |w<22|)- 
One can see from Eq. ( p8| ) that the latter can be obtained only at the price of a cancelation between 
x\ and y 2 , and/or x 2 _ and y 2 _ ; in this case v\ and v 2 form a pseudo-Dirac pair. 

As far as \m\2\ = x\X- + yiy- is concerned, one can use the fact that, generically, i.e. barring a 
close correlation between the z's and y's, one has |y + | ~ ~ max ( 1 2/2 1 ; 1 2/3 1 ) , and analogous relations 
for the x components. Now from ( |36|) one gets the bound 

|mi 2 | ~ \x\x+ + yiy + | < |LT e3 |m@ (39) 



that can be used in (34) to test the large mixing mechanism in hierarchical models as discussed below 
in more detail. 

If m@ is dominated by one right-handed neutrino and the solar mixing angle is due to some 



cancelation between the two terms in the denominator of Eq. (31), there are two situations as far as 
the two remaining neutrino mass eigenvalues are concerned. They can exhibit some hierarchy, with 
m 2 — m Q 3> m i, or they can be comparable in magnitude (but not much larger than ?n@ to keep 
large mixing angles stable under the RGE evolution). Let us discuss the first possibility in detail and 
then comment on the second one. 
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3.1.1 Double Dominance 



In the presence of the large solar angle, the most natural realization of a hierarchy between 771,2 and mi 
involves a dominance mechanism analogous to the one assumed in the atmospheric sector. Indeed, if 
xi *C yi, X- <C y_, the solar neutrino scale is dominated by the right-handed neutrino with mass M y , 
namely mg ~ 777,2 ~ \yf + y 2 | (the choice of x negligible with respect to y is not restrictive: yi <ii, 
y- <C X- corresponds just to the exchange M x <-> My). When xi,x- are much smaller than yi,y~ 
the formulae above simplify considerably also because sines and cosines factorize out. Indeed, it turns 
out from the condition ([}(]) that 02 — X ~ — ( a i + x) = 0- The CP violating phase x is 

X ~ arg(y_) - arg(yi). (40) 

The expressions for the solar angle and the solar mass scale become simply 

tan#i2 « r , m « |y_| 2 (l + r 2 ) , r = j^-y . (41) 

Then, as was the case for the atmospheric angle, a large solar angle results if the ratio of two couplings, 
yi/y-, is of order one. So, the prediction for a large solar angle is stable in the sense specified in the 
introduction. Together with the assumed dominance in m@, this corresponds to a double-dominance 
pattern (or triple-dominance pattern, since the lightest mass eigenvalue, mi, is then dominated by 
the third right-handed neutrino). LMA requires 0.5 < r < 0.9 and |y_| 2 = (3 - 8) x 10~ 3 eV, while 
LOW requires 0.7 < r < 1.1 and |y~| 2 = (1 - 3) x 10~ 4 eV (using the 99% C.L. intervals of Ref. gj|). 
A maximal solar angle would require a tuning of r to 1. Finally the Major ana phase in V associated 
with m-2 is A2 = arg(V22) ~ arg(yi). In the limit mi — > 0, there is no other CP-violating phase; in 
general however U contains a third CP-violating phase associated with mi, Ai = arg(Vn) (in our 
phase conventions, arg(V33) = 0). 

A very constraining conditions holds for the double dominance pattern which follows from the 
©, © and © 

\U e3 \ > I s in2# 12 ^ . (42) 
2 m@ 

For LMA this means that U e s should be within one order of magnitude from its present experimental 
limit, while it can be much smaller for LOW, as shown in Figure ^ a). 

We finally note, from the assumed perturbativity property of the Yukawa couplings, that M z < 
5 x 10 14 GeV, while M y < ^| x 5 x 10 14 GeV. For the LOW solution M y can lie above the unification 
scale (when the corresponding Yukawa couplings are of order one). 

The double dominance pattern can also accommodate a small solar angle, with yi <C y_ instead 
of yi w y~. For the (now strongly disfavoured by the data) SMA solution, one needs 0.01 < r < 0.03 
and \y~\ 2 = (2 — 3) x 10 -3 eV. The CHOOZ angle is constrained to be small; more precisely, it lies in 
the following range: 

tanfl 12 ^ < \U e3 \ < tan0 12 . (43) 

m@ 

Thus |?7e3| is at least one order of magnitude smaller than its present experimental limit; this is to be 
contrasted with the large angle solutions for which there is no theoretical upper bound on \U e3 \. 
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a) Double dominance 



b) Pseudo Dirac 
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Figure 2: Lower limits on U e s for double dominance a) and pseudo-Dirac b). To sketch the regions of 
the LMA and LOW solutions we have shown in blue the 95% C.L. contours of the first reference of 
0. The lower (upper) curve displayed for each color corresponds to nig = 2 (4) x 10 -3 eV 2 . 

3.1.2 Pseudo-Dirac Solar Neutrinos 



The alternative option for generating a large solar angle is to assume a pseudo-Dirac structure with 
\mn\, I m22 1 < I m i2|, i-e. 



I —2 i — 2 1 I 2 i 2 

1 2/1 , \y_ +x_ 



•c \yiy- + x\x-\ = m. 



(44) 



This pattern leads to a strongly mixed pseudo-Dirac pair of neutrinos with masses 772-2 — —mi ~ 777; 
\U e 3\m@ > fh 3> 77i0. As already discussed in Section ^T|, the conditions fl44|) require cancelations 
between the y's and the x's, which is most easily implemented with a pair of pseudo-Dirac right- 
handed neutrinos, with M y ~ — M x ~ M, together with some hierarchy among the Yukawa couplings. 

Since a large angle is an automatic consequence of the smallness of |mn|, 1 777,22! implied by 



no particular correlation between them is required, and the ratios in the test relation (34) are 0(1). 
Replacing |mi2| by its upper bound |£7 e 3|m@ one finds the constraint 



m@ tan 2 6 



12 



^5 \Ue2,\ 

1 — tarn 6/12 



(45) 



As shown in fig. § b), in the LMA region \U e ^\ should be very close to its experimental bound and 
saturate the lower bound in Eq. There is more room for the LOW solution. The lower bound 



on \U e ^\ in Eq. (|4q) depends on how #12 is close to tt/4; the smaller 1 — tan #12, the larger \U ( 



t'3 
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From the perturbativity property of the Yukawa couplings one can derive an upper bound on the 
right-handed pseudo-Dirac neutrino mass, M < x 5 x 10 14 GeV. Since > | t/ e 3 1 1 , the pseudo- 
Dirac right-handed neutrino pair could be heavier than the third right-handed neutrino, whose mass 
satisfies M z < 5 x 10 14 GeV. 



3.1.3 Two-scale spectrum: rri\ ~ ra^ 

The double dominance mechanism gives a fully hierarchical mass spectrum, m\ <C mi <C mz - However, 
m\ ~ m,2 is in principle a viable option too. This corresponds to the case where no single right-handed 
neutrino dominates in the left upper 2x2 submatrix in Eq. (|38|). Then, full Eqs. ( |30|) to (|32^ ) must be 
used, and there is no simple correspondence between the solar neutrino parameters and the couplings 
of one right-handed neutrino as in the case of double dominance. The large solar angle requires 
\rri22\ cos(a2 — x) ~ \ m n\ cos(ai + \) an d \?ni2\ ~ |"i<ii| 5 |"i22|, (which is obtained for (i) x\ < y~\ ~ y_ 
if x_ < y_; (ii) yi < x\ ~ x_ if y_ < x_; (hi) xi < y x ~ y_ or y x < x\ ~ x_ if x_ ~ and the 
solar neutrino scale is approximately given by m & ~ \y^_ + x 2 _\ cos(a2 — x) + + ^il cos(ai + x)- 
Thus, sensible variations in the magnitude of 612 and m are expected by slightly perturbing the many 
parameters entering the previous expressions, in particular the phases. Indeed, as already mentioned, 
the phases work against a large solar angle (they are not dangerous in the double dominance case 
because the cosinus in the denominator of eq. (^) factorize out). Thus, the two-scale spectrum 
cannot arise from the see-saw in an economical and robust way. The CHOOZ angle satisfies again the 
same lower bound as in the double dominance case: |t7 e 3| > O {m & /m®) . However, due to the large 
number of parameters involved, this bound is less strict than Eq. (f42|). For completeness, note that 
the SMA solution would be accommodated with (i) y\ <C y_ and x\ < y_ if x_ C (ii) x~i <C x_ 
and yi < x„ if y_ <C x_; (hi) x\ <C x_ and y\ <§C y_ if x_ ~ y_. The solar neutrino scale is then 
given by rriQ ~ \y^_ + x 2 _\ 2 — \y\ + x\ | 2 , and the CHOOZ angle lies in the same range as in the double 
dominance case if both x\,y\ -C max(x_, y_). 



3.2 Inverted Hierarchical Spectrum 

Let us now discuss the alternative case of an inverse hierarchy in the neutrino mass spectrum. We 
shall work with the variables (xi , ui , v-i), introduced in section ( |2.3| ) such that L is given by (]24|). Up 
to irrelevant phase redefinitions, the implementation of the atmospheric neutrino oscillations require 
v\U-\- ~ m@ , with the other components much smaller. As before, the unphysical phases W can be 
chosen so that #23 and #13 are real. Expanding in the smaller matrix elements of L, one obtains the 
following expressions for these angles: 

tan 9 23 « - ( =± ) + ~ V - + X1X ~ , (46) 







U3 _ 1 




U 2 ' 










(vx 2 


+ 




Vl 





tan 13 « -- + ( ^ ) U ~ V+ + X ~ X+ « - V - 1 + ™± )+ X -^± (47) 
where the variables (x±, v±, u + ) are defined in (|l5|), with 8q = arctan(— U3/1L2), U- = 0, and 



xi^i = (vix\ + v _x_) , x_v\ = (uix_ — v_x\) , u\V\ = , 

u_vi = —V-Ui , V- = , v\ = v\ + -y 2 , n + wi « m@ . (48) 
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In this approximation, the matrix m defined in (|29| ) reads 

(2u±vx + x\ v\u+ + u\v + + 

t>i7i+ + + 2u + u + + x 2 + | . (49) 

x 2 _ 

By replacing these matrix elements in (|30|) , (|3l|) , ([32^) and (|33|), one obtains the expressions for tan 2 #i2 
and itiq. 



boost given in (p5|), it is worth defining the rapidity A = \ log(u + /t>i), such that 



Since the matrix m and the oscillation parameters discussed in this section are invariant under the 

1 
2 

A /=— „. _ „-A 



u + ~ e ui ~ e y/m@ . (50) 

while the parameters (ui/u+), /fi), (v-/vi), are A-independent. They appear in the approximate 
expressions: 



2m@ u + ui m@ 

9 „ U\ v + x 2 — x\ , , 

tan 2 012-1 ^ 2— - 2— + — ± , (51) 

u + t>i m@ 

where the phase dependence discussed in Section ^ has been omitted for simplicity. 

The important feature here is the relation |mi2| ~ v\u + ~ m@. Since in this inverse hierarchy 
scenario, the large solar neutrino mixing is related to the pseudo-Dirac structure, we preclude a strong 
cancelation in the right-hand side of the relation ([54|) , namely a precise cancelation between the matrix 
elements m\\ and 77722 of (|49]). In particular, the two expressions in ( |5l|) must be of the same order. 
Thus, we derive from (|34|), by replacing 1 77112] by m@, the constraint 



,2 



;i-tan 4 0i 2 ) 



mk tan 2 12 



~ ^ T . (52) 



The comparison with the fits to neutrino experiments is displayed in Fig. ^. Therefore the inverse 
hierarchy patter is strongly disfavoured for LMA unless we complement the pseudo-Dirac pattern with 
a fine tuning of m\\ and 77122- As for the LOW solution one needs the solar mixing angle to be very 
close to its maximal value, namely, |1 — tan 2 #12! < 10 -4 . 

Generically, i.e. without fortuitous cancelations, one expects that these strong upper limits on m\\ 
and 77722 will translate into correspondingly small values for |v + /ui|, |ui/u_|_| and \x\ |/m@, |x^_|/?77@ . 
As already argued, we assume the natural relations ~ |«+| and ~ |^+|, to obtain a limit 

lOesI < |4 (53) 

which is below the present experimental limits for LMA and well below it for LOW. Notice that this 
conclusion is quite independent of the tuning to get m & discussed above. 

Actually, one can distinguish two typical mechanisms to implement the two small scales tuq and 
7773 ~ x 2 as well as the small deviation from maximal solar mixing angle #12: 
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Figure 3: The curve determined by Eq. (52). The blue regions are like those of Fig 



1) Double dominance by a pseudo-Dirac pair, so that mg > 772.3: when the contributions of the terms 



of the form be neglected in (|4q), (|47j ) and (|5l[). Both the atmospheric and the solar neutrino 

oscillations are basically controlled by the sector (u , v), hence tuq and 612 are dominated by the 
right-handed pseudo-Dirac pair contributions. In this case, the ratios {u\/u + ), (v + /v\), (v-/vi), are 
approximately determined by and (|5l|). 



2) Single dominance of vtiq , with ~ 7713: when the contributions of the terms of the form mvj can 
be neglected in (f46D, ( ^7| ) and ( |49| ) (but for u + v±), so that the terms XiXj control the solar parameters 
mo and the deviation of 612 from 7r/4, with a dominance by the third right-handed neutrino. Here, 
Xr, x — , X-\- , are approximately fixed by (^7|) and ([H|). 

As we now turn to discuss, charged lepton flavour violating (CLFV) decays can provide some 
discrimination between these different possible patterns. 



4 Charged lepton flavour violating radiative decays 

In the previous section, we identified the patterns for L that reproduce in the most natural way, from 
a bottom-up point of view, the experimental results. Of course, neither the absolute scale nor the 
hierarchy of right-handed neutrino masses can be deduced from oscillation data. Other experiments 
are needed to gain further insights; in particular, searches for flavour-violating charged lepton decays 
such as r — > /J7 and /i — ► ej offer the possibility to test another combination of Y and M than the 
one associated with oscillations. 
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It is well known that the flavour structure of soft terms in the slepton sector can lead to strong 
violations of lepton flavour in supersymmetric extensions of the Standard Model. Even if one assumes 
that the mechanism of supersymmetry breaking is flavour-blind, with m 2 ~ . . = ?rt|- • = m^Sij and 
Afj = AoYfj at some scale Mjj ~ Mpi, the soft terms can receive additional flavour-dependent contri- 
butions from various sources, which may lead to an observable rate for processes such as r — > fij and 
fx — > erf. Among the possible contributions are the radiative corrections induced by the right-handed 
neutrino couplings Y^i [^] , radiative corrections in the context of grand unification^] and other possible 
contributions from unknown physics between Mqut and Mp\. 




Figure 4: Upper bound on - — from the present experimental limit BR(t — > /^y) < 1.1 x 10 6 as 

a function of the mean sneutrino mass m,y and of the SU{2)l gaugino soft mass M2, for tan/3 = 10, 
Aq = and sign(//) = +. The numbers in brackets correspond to an improvement by a factor 10 3 in 
the experimental limit on BR(t — > /ry). The dark gray region is excluded because the lighter stau 
would here be the LSP, while the light grey region corresponds to ttiq < 0. The region below the 
dash-dotted line could be explored by searching for LFV slepton decays at future colliders p8j . 

The running of the slepton masses from Mjj to the scale where right-handed neutrinos decouple, 
due to loop corrections involving the Yukawa coupling matrix Y, induces flavour non-diagonal terms 

4 For instance it is well known that in 5(7(5) an additional source of LFV are the Yukawa couplings of the colored 
triplet @. 



18 



1000 



800 



600 



400 



200 




200 400 600 
M 2 (GeV) 



800 



1000 



Figure 5: Upper bound on 



from the present experimental limit BR([i — > cy) < 1.2 x 10 as a 



function of mo and M2, for tan/3 = 10, Aq = and sign(/x) = +. The numbers in brackets correspond 
to an improvement by a factor 10 3 in the experimental limit on BR(/i — ► e~f). The regions in dark 
grey and light grey are defined as in fig. 1. 



in the slepton mass matrices which are proportional to the off-diagonal entries of 

C = y f In (Mu/M) Y . 



(54) 



These flavour-dependent soft terms induce LFV charged leptons radiative decays l{ — > lj*y via loops 
of charginos/neutralinos and sleptons. The corresponding amplitudes have been fully computed, in 
the mass insertion approximation, in Ref. [39]. For moderate and large values of tan/3, the dominant 
contribution arises from loops where charginos and sneutrinos circulate, with an insertion of the off- 
diagonal element of the sneutrino mass matrix m| ^ ; . This yields a branching ratio 



BR(k — Z i7 ) 



a: 



G 2 f{M 2 ,H,mo) \m\..\ 2 tan 2 13 , 



(55) 



where / is a function of the SU(2)l gaugino mass parameter M2, the supersymmetric Higgs mass \i 
and the mean sneutrino mass m^. Thus the dependence of BR(^ — ► L7) on the seesaw parameters 
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is encoded in the off-diagonal entries of the matrix C; this remains true when the subdominant 
contributions are taken into account. 



In figs. U and ^j, we show the upper limits on (m~) ^ r / m~ and (m 2 ) eiM /m 2 that can be inferred 
from the present bounds on the branching ratios for r — ► fij and \i — > cy, respectively. The bounds 
are displayed in the (m/^-Ma) plane, the variables which they are mostly sensitive to, for tan/3 = 10, 
Aq = and sign (/x) = +. However, since BR(Zj — > Ijj) scales as tan 2 f3 for moderate and large values 
of tan (3, the upper limits associated with other values of tan (3 can be obtained upon multiplication 
by (10/ tan/3). We also give in brackets the limits corresponding to an improvement by three orders 
of magnitude in the sensitivity to the branching ratios. Such an improvement is indeed expected for 
fi — ► e7 |4C]. Prospects for r — > lij are currently less optimistic, but searches for the LFV decay 



X2 ~^ XiA 47 " a t future colliders could provide limits on (m~) « T / m~ of that order of magnitude in the 
region of the (ms, M2) plane indicated on figs. |I| and || p8fl . We do not consider the process r — ► e7, 
which is much less constrained than /j — > e7 by experiment, and hence does not provide additional 
information on the seesaw parameters. 

In the framework of mSUGRA the coefficients are related to the (m-^r/m- by 

2 „ _ 3m| + _A| 

L 2J 8vr 2 J v 7 



The limits in figs.|| and corresponds to those in figs. |] and ||. The latter are taken from \21\ where 
a more complete discussed can be found. 

In view of the following discussion, it is useful to express the quantity Cji in the parametrization 

Cjt h ( ^ + ^ ln ( f) . (57) 

At this stage, nothing has been assumed about the structure of L, and Eq. ( |57[ ) is completely general. 
We shall neglect CP violation in the following, and therefore replace ZjZi by ZjZi, etc, with an additional 
— sign if the corresponding right-handed neutrino has a negative CP eigenvalue. 



4.1 Hierarchical mass spectrum 

In the case of a hierarchical light neutrino mass spectrum, m@ ~ z 2 ., and it is convenient to rewrite 
Eq. (H|) as 

z*Zi M z (M V \ y*Vi M y (Mjj\ x* Xi M x ( M v \ 

where M@ = v 2 /m@ ~ 5 x 10 GeV. The dominance hypothesis, together with the assumed per- 
turbativity of Yukawa couplings, implies M z < M@. In the case of double dominance, one also has 
My < M Q = v 2 /m Q (with typical values M ~ 5 x 10 15 GeV for LMA, and M ~ 10 17 GeV for 
LOW), and M x < Mi = v 2 /m\ (M\ S> M ). Since mi is not known, the constraint on M x is not 
very informative; however if one gives up the dominance in the solar neutrino sector, one has both 
M x < Mq and M y < M Q . In the pseudo-Dirac case, Eq. (p7|) further simplifies to 

^ Mn_ (Mu\ yjVi - xjXi M_ (Mu\ 
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Figure 6: Upper bound on C23 = (Y^ln(-^)Y)23 from the present experimental limit BR(t — > /ry) < 
1.1 x 10~ 6 as a function of the mean sneutrino mass m^ and of the SU(2)l gaugino soft mass M2, 
for tan/3 = 10, Aq = and sign(/x) = +. The numbers in brackets correspond to an improvement by 
a factor 10 3 in the experimental limit on BR{r — > /17). The dark gray region is excluded because the 
lighter stau would here be the LSP, while the light grey region corresponds to m\ < 0. The region 
below the dash-dotted line could be explored by searching for LFV slepton decays at future colliders 

Note that the \n.(Mu / M) term in Eq. ( |59|) is proportional to y^yi — xjx^, not to yjyi + XjXi. This 
is important since cancellations may occur in one of the two combinations, depending on how the 
pseudo-Dirac structure is realized. More precisely, one may have \y1y2 — x±X2\ <C \y1y2 + ^1^2! if 
cancellations occur both in mu and 77222 (implying a suppression of the contribution of the pseudo- 
Dirac right-handed neutrino pair to C12), and \y2y3 + ^2^3) "C \y2y3 — X 2 X 3\ if cancellations occur in 
77722- Furthermore one has M < v 2 /rh ~ y^2|1 — tan 2 9\2\Mq. Note however that the actual upper 
bound on M can be smaller, even in the presence of Yukawa couplings of order one. In fact M ~ 
a/211 — tan 2 <9 12 | M Q requires, in addition to Yukawa couplings of order one, both cancellations in mu 
and m 2 2, and \y\\ ~ \y~\). 
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Figure 7: Upper bound on C\i = (Y^ln( jf-)Y)\2 from the present experimental limit BR(/j, — > ej) < 
1.2 x 10~ n as a function of and M2, for tan/3 = 10, Aq = and sign(/x) = +. The numbers in 
brackets correspond to an improvement by a factor 10 3 in the experimental limit on BR(fi — ► e7). 
The regions in dark grey and light grey are defined as in fig. 1. 



4.1.1 r — ► yrf 

In the double dominance case, the coefficient associated with the decay r — > /ry is 

1 ■ no M * 1 ( MjJ \ V*V* M V 1 ( Mu \ X2X 3 M * 1 ( MjJ 



C 23 ^ - sin20 23 In ^ i^^lnp + In ^ . (60) 

2 M a v M *y 4 M @ \ M y) 4 M @ V M */ 

Although, from Eq. (f4l|), 1 2/22/3 1 ^5 m © an d 1^2X3 1 <C m , the second and third terms in Eq. (60) can 
be as large as ln(M[//M ) and ln(Mrj/Mi), respectively, if M y and M x are close to their respective 
upper bounds. In the pseudo-Dirac case: 

1 . M z (Mu\ 2/22/3 ~ x 2 x 3 M (M v \ 

C 23 - -sm20 23 — In i^—j + ^ — In . (61) 

In general one expects the contribution of the pseudo-Dirac pair of right-handed neutrinos to be small, 
since from Eq. (f£3|) |2/22/3+ 2; 2£3| <C m / y/2\l — tan 2 6*12 1 - Still it can be of order one when cancellations 
occur in 777-22 due to both 7/3 ~ ±1x3 and 7/2 — ±i&2, since in this case 1 2/22/3 — %2%3\ 3> 1 2/22/3 + ^2^3!- 
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If M z is the heaviest right-handed neutrino mass, Eqs. (|60| ) and (|H|) reduce to C23 ^ \ sin 2#23 
\n(Mjj/M z ); in this case the observation of the process r — > fij would amount to a measurement of 
M z as a function of M2 and m v - assuming no other contribution to the LFV slepton masses, e.g. 
from the supersymmetry breaking mechanism itself. In general, however, M z can be smaller than 
M x or My (resp. M), and the second or third term may give the dominant contribution to C23. 
Hence a measurement of BR(r — > /17) can only provide an upper bound on M z . A small branching 
ratio, I C23 1 <C sin 2#23 hi(M(//M@), would point towards two alternative scenarios |2l[]: (i) a scenario 
characterized by small Yukawa couplings in the neutrino sector and a "low" lepton number breaking 
scale, with both M z <C M@ and M x ,M y <C M (resp. M <C y/2\l - tan 2 6 12 \ M ); (ii) a scenario 
characterized by the dominance of a lighter right-handed neutrino in the atmospheric sector and a 
"high" lepton number breaking scale, with M z <C M@ and M y ~ M or M x > M (resp. M close to 
\/2\l — tan 2 #12! M ), the contribution of M y (resp. M) to C23 being suppressed by IX2X3I, 1 2/2J/3 1 <C 
m (resp. 1 2/22/3 — £2^3! <C m /y / 2|l — tan 2 6*12!). In the latter, the B — L symmetry is broken close to 
- or even above - the scale where the MSSM gauge couplings unify (i.e. 2 x 10 16 GeV), which favours 
grand unified theories based on SO (10) and larger gauge groups. On the contrary the former appears 
to disfavour such theories. 



4.1.2 fi —> e7 

Let us first discuss the LMA and LOW solutions, which are strongly favoured over SMA by the solar 
neutrino data. In the double dominance case, the coefficient associated with the decay [i — > ej is 

— sin0- — i« (— ^ _l ll^I Ml i„ ( Mu \ jl ^1 Ml i„ ( Mu 
In the pseudo-Dirac case 



C ~ - s y 23 -i n p | ELLL± m_ l n l^i + ZiZ£ ^ hi . ( 62 ) 



C 12 - * sin.23 ^ In f ^ + ^ 2 -^ 2 # m f ^ ] . (63) 



One can show that the contribution of the pseudo-Dirac pair of right-handed neutrinos is always sup- 



pressed relative to ln (Mjj /M Q i/2|l — tan 2 9i 2 \ ), due to either 1 2/12/2 — ^1^2! "C ^o/v^P- ~~ * an2 ^12 1 
orM< \/2|l — tan 2 6*12 1 M , while in the double dominance case M y or M x could in principle give a 
contribution of order one to C\ 2 . 

Expectations for BR(,u — > e7) are very model-dependent. It is useful for the discussion to distin- 
guish between two regimes differing by the value of \U e 3\. In the "large |C^ e 3|" regime, characterized 
by 2|i7 e 3| ^> sin 2#i2 m /m@ in the double dominance case (resp. \U e ^\ 3> m & / \J2\1 — tan 2 6\ 2 \ m@ 
in the pseudo-Dirac case, where this regime exists only for the LOW solution), the first term in Eqs. 
(|62^) and (|63|) yields a Mj-dependent lower bound on | C12 1 , 

|C d >sin« 23 tannin (^|) . (64) 

Eq. (|64|) can be used to convert an experimental limit on BR(/i — ► e^) into an upper bound on M z as a 
function of M 2 and m v - assuming \U e 2,\ is known. Actually the present limit on BR(/i — > e^f) already 
provides strong constraints on this regime; in particular a large M z , M z ~ M@, is already excluded 
over a large portion of the (m u , M 2 ) plane, especially for LMA. The contribution of the other right- 
handed neutrinos also yields non-trivial constraints. In the double dominance case, | C 12 1 <C 1 requires 
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1 2/12/2 1 sin20 1 2m Q /2 and/or M y <C M Q (there are no analogous constraints for IX1X2I and M x since 
mi is unknown). In the pseudo-Dirac case, \C\2\ "C 1 requires \y1y2 — X\X2\ <C r^/y^ll — tan 2 12 | 
and/or M <C \/2|l — tan 2 6>i2 1 M & , which is easily obtained due to the pseudo-Dirac structure. With 
the expected improvement by three orders of magnitude in the experimental limit on BR(/j — ► ej), 
those constraints should become rather stringent - unless /i — ► e^ is observed. 

Note finally that if atmospheric neutrino oscillations are dominated by the heaviest right-handed 
neutrino mass (i.e. M x ,M y < M z , resp. M < M z ), the contribution of M z dominates both in 
BR(r — > and in BR(/i — ► e7). This leads to the prediction 

|C^| „ tang 13 

|C 23 | ~ COS023 ' 1 ' 

While Eq. (pq ) remains valid over a large region of the seesaw parameter space with M x > M 2 or 
My > M 2 (resp. M > M z ), a deviation from this relation is possible only if M z is not the largest 
right-handed neutrino mass. In this case one may have \C\2\ 3> \U e z\ m (Mu/M@) ~ 8 \U&\ (which is 
not excluded for the LOW solution yet), hence an observable rate for /j — > erf. 

The "small |C/e3|" regime, characterized by 2|C/ e 3| ~ sin20i 2 m0/m@ in the double dominance case 
(resp. \U e -^\ ~ m Q /y / 2|l — tan 2 6*12 1 m@ in the pseudo-Dirac case), differs from the large \U e "A regime 
by the absence of a lower bound on | C12 1 : here z\ — > is perfectly compatible with the existing data. 
Thus BR(/i — ► ej) can be extremelly small even though M z ~ M@ (i.e. even though BR(r — > fij) is 
large), unlike in the previous regime. However | C12 1 3> | t^e3 1 ln(Mjj/M@) is also a possibility when M x 
or My S> M@ > M z (resp. M S> M@ > M z ), and /1 — > e^ should be observed at PSI in this case. 

The discussion of the (strongly disfavoured) SMA solution follows the same lines. Again one can 
distinguish between two extreme regimes. In the "large |J7 e 3|" regime, where | C/e3 1 ~ tan #12, the 
lower bound (^4|) holds like for the large angle solutions, as well as Eq. ( |65[ ) in a large portion of 
the seesaw parameter space including the case where M z is the largest right-handed neutrino mass. 
Note however that BR(/i — ► cy) cannot be as large as in the LMA or LOW case, due to the upper 
bound |Ci 2 | < tan #12 ln(M[//M@) (to be compared with |C 12 | < ln(Mu/M Q ) for LMA/LOW). In the 
"small I C/" e 3 1 " regime, where \U e s\ ~ tan#i 2 m0/m@, BR(/x — ► ej) can be extremelly small without 
conflicting with existing data. Again the observation of fi — > e-y with |Ci 2 | 3> \U e $\ ln(M[//M@) is 
possible if one of the two lightest right-handed neutrinos dominates in the atmospheric sector, i.e. if 
M x or M v > M@ > M z . 



4.2 Inverted Hierarchy Mass Spectrum 



When the neutrino mass matrix is endowed with an approximate pseudo-Dirac pattern as required by 



an inverted hierarchy in its eigenvalues, it is convenient to rewrite the coefficients Cjj given by (58) in 



terms of the variables defined in Section 3.2 as follows: 
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(66) 



Since, as discussed, an approximate pseudo-Dirac texture in the neutrino mass matrix is more naturally 
realized by a corresponding pseudo-Dirac pair of right-handed neutrinos, we have associated the masses 
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±M to this pair. The dependence on the rapidity A shows the dependence on the boosts (|25| ) and 
displays an example of the dependence of the Cj,- on the orthogonal tranformations discussed in Section 




The Yukawa coupling associated to the largest entry in the matrix L, which in this case must be 
either u+ or v\, is perturbative if M < e _2 ' A 'M@ . The best limit on M x is obtained when the xi 
contributions saturate the quantities given by ([17]) and d5l|): M x < {m^/2U^m^)M@ < M@/\U e s\ , 
where the last bound relies on (|53|), so that M x lies well above Mq. 



4.2.1 t — > /17 

The coefficient associated to the r — ► ^7 decay is 

1 9> M {Mu\ x\ + x 2 , ,M X ( Mu\ , . 

where tan 7 = Since sin 2023 is large, the first term gives a contribution 0(1) when M 

approaches its bound e _2 ' A 'M(Q , which could be much below the gauge coupling unification scale if 
A > 1 . On the contrary, this contribution is suppressed by a factor e~ 2 ' A ' for A < . As for the second 
term, it is strongly suppressed by the bounds on the factors, such that (x 2 + + x 2 _)/2m@ is expected 
to be 0(m|/2m|). This suppression can be compensated if the mass M x gets close to the bounds 
imposed by the Yukawa coupling perturbativity condition. For LMA the bound on M x lies near the 
gauge coupling unification scale, but for the more appropriate LOW solution it is around the Planck 
scale, which is clearly nonsense. Therefore, some suppression is expected on quite general grounds for 
this term. 



4.2.2 fi -> e7 



The coefficient C12 corresponding to the fi — > e^ decay is 
C\2 ~ 



x±\ x 



+ 



^ e 2X + v_± e -2A \ S 23 + v - e -^ sin e 23 

U + Vi ) Vl 



x - . (a , s M x (M V 
sm(&23 + 7)^^ In 



Mi 



(68) 



For simplicity, let us separately consider the two typical situations previously identified in Section 



1) Double dominance by a pseudo-Dirac pair; replacing the ratios (ui/u+), (v+/vi), (v-/vi), by the 
expressions obtained from (|47| ) and (|5l|), yields 
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4m 2 



Y cosh 2A + (tan 612 — 1) sinh 2A ) cos #23 — e 2A tan #13 sin 9 



703 



M (Mv 
M 



Mm 



O 
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(69) 
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where the last evaluations corresponds to the generic case - i.e., to a natural pattern in the sense 
adopted in this paper - and A > . This limit is very small for the LOW case, the one most relevant 
here. Notice that in this case, the approximate relation holds: C12/C23 ~ O (m0/2m|) The upper 
limit on C12 can be saturated only if M is close to its bound which implies a relatively large r — ► fij 
decay. Of course, this contribution can be overcome by the other one if M x is very large. 



2) Single dominance of m© ; then the components Xi can be evaluated from ( f47j ) and ( |51[ ) in terms of 
777-0, (tan 2 612 — 1), £7 e 3. However, the generic relation: x\x^ ~ x\x + < [x\ + x\)/2 m m Q /4m@ leads 
to the estimate 

^ ~ o (&) £ m ■ (70) 

\4mlJ M @ \M X ) v ' 

Therefore, in this case, C12 could be of 0(1) only if M x lies close to its bound M@/\U e ?\ . However, 
for the more suitable LOW solution, ( |53| ) implies that M x is too close (or above) Mjj for the result 
to be sensible. By assuming that the contributions from the x^ dominate in both C\i and C23 their 

ratio is given by: C12/C23 ~ x\j \J x\ + x 2 _ cos(#23 + 7), which is less suppressed than in the previous 
case and just tells us that the right-handed neutrino with mass M x is more coupled to e than to \i 
or r. This is not possible in the previous case of double dominance by a pseudo-Dirac pair because 
this ratio is suppressed by the conditions to implement a pseudo-Dirac pair. Of course the important 
prediction from the experimental viewpoint is for each coefficient. 



5 On the Physical Interpretation of R 

The ambiguities in the extraction of the see-saw parameters from the measurements of M. v can be 
expressed in terms of the complex orthogonal matrix R defined by: 



R I ' I U ] = M~ x l 2 Yv = L (71) 

\frh~i, 

The matrices are defined as in the previous sections. The nice property of this choice of R is its 
immediate interpretation as the dominance matrix in the sense that: 

(i) R is an orthogonal tranformation from the basis of the left-handed leptons mass eigenstates to the 
one of the right-handed neutrino mass eigenstates; 

(ii) if and only if an eigenvalue m, of A4 U is dominated - in the sense already given before - by one 
right-handed neutrino eigenstate Nj, then \Rji\ ~ 1 ; 

(Hi) if a light pseudo-Dirac pair is dominated by a heavy pseudo-Dirac pair, then the corresponding 
2x2 sector in R is a boost. 

It is instructive to illustrate this property by a simple example, which in a sense also motivates 
the double dominance assumption in Section [||. Consider a hierarchical neutrino mass spectrum, 
TO3 ~ m@ ^> ni2 ~ nriQ m\ . For the sake of the example, we fix the other 9 parameters (all 
phases are neglected) as follows: (a) The Y eigenvalues are very hierachical with 1/3 > 1/2 > 1/1 ; 
(b) bimaximal mixing angles in U; (c) the Vl matrix defined, as usual, by the diagonalization Y = 
Vjjdiag(yi, 1/2, Vz)Vl, is very small or very close to the deviations of U from the bimaximal case. This 
last assumption is interesting because it selects the case where the charged lepton flavour violations 
are very small, while the mixings are maximal in M u . Then we get the approximate expressions for 



2(3 



the M eigenvalues: 



M, 
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Mo 



M, 
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4mi 



(72) 



where M3 is much larger than M2 <C M@, while Mi <C M Q depends on the hierarchy in the Y 
eigenvalues. Notice that our assumption of hierarchical spectrum for the left-handed neutrinos and 
small charged lepton flavour violations leaded to a pattern of M eigenvalues of the kind discussed in 

m. 

As for the dominance matrix R, which in this case is a rotation matrix, we get 



R 





(73) 



The entries (close to) 1 in R clearly indicate that m@ is dominated by N2, m & by N\ and, of course, 



m\ by A3 . In comparing with the discussion in Ref . [ 21 1 , one should take into account the different 
ordering of the right-handed neutrinos (by increasing masses) therein. It is worth stressing that our 
assumptions of hierarchical light neutrinos and Y eigenvalues yielded the double dominance pattern 



in (72). Analogous exercises show how the dominance of a pseudo-Dirac pair displays in R. 



As a second, more familiar example, we consider the class of models already mentioned, in which 
the fermion mass matrices are explained by an abelian flavour symmetry with only positive charges 
for all leptons. It is well known that in this case A4 U is independent of the right-handed neutrino 
charges and the magnitude of its matrix elements {Ai u )ij are of 0{e h+l ->). It is easily seen that, in 
this framework, the R matrix elements come out all of 0(1) , characterizing the fact that all heavy 
neutrino states are comparably contributing to all light neutrino mass eigenvalues. 



6 Outlook and Conclusions 

To have a glimpse at the large solar mixing angle from a bottom-up perspective, we adopt the re- 
quirements of robustness and economy and we study the patterns for Y and M which, in this sense, 
naturally account for neutrino phenomenology. A quantitative appraisal of the structure of subleading 
terms and of the effect of CP violating phases, is done by means of reliable approximate expressions. 
Three scenarios, which admit a suggestive physical interpretation in terms of the role played by right- 
handed neutrinos, are favoured: single dominance in the atmospheric sector accompanied by single 
or pseudo-Dirac dominance in the solar sector, in the case of a hierarchical neutrino mass spectrum; 
pseudo-Dirac dominance in the atmospheric sector in the case of an inverted hierarchical spectrum. 

The phenomenology of the three scenarios is quite rich. In the hierarchical case, experimentally 
relevant lower bounds for U e z are derived. In the inverted hierarchical case, on the contrary, upper 
bounds of the order of permille are obtained. All three scenarios predict a rate for (3f3Qv too small 
to be observed (for a recent exhaustive analysis see Ref. j4l|| ). Connections with the rate of lepton 
flavour violating decays r — ► fij and /x — > e^f have also been discussed. These rates are linked to the 
nine hidden parameters of the see-saw. The strong dependence on the three right-handed eigenvalues 
is manifest, but also the parameters belonging to the R matrix play an important role. For instance, 
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in the pseudo-Dirac dominance scenario, there is a very strong dependence on the 'rapidity' defined 
in Eq. (^5|), which represents the magnitude of one of the boosts contained in R. 

It is worth emphasizing that, while being 'minimalist' from the bottom-up perspective, these 
favourite scenarios are inconsistent with a simple flavour symmetry like a U(l) with only positive 
charges for leptons. Indeed, they suggest the presence of much richer flavour symmetries, like £/(l)'s 
with holomorphic zeros from supersymmetry [13, 14] or non-abelian flavour symmetries [15]. Finally, 



we also stress that the matrix R [18] admits a suggestive physical interpretation. Namely, it encodes 
the informations on the dominance mechanism at work. 



Note added in proof: As stressed in the Introduction, the different dominance scenarios analysed 
in this paper from a bottom-up perspective have been discussed several times in the literature, mostly 
from a top-down viewpoint. The dominance of one right-handed neutrino first appeared in Ref. |25| 1 
and both the single and double dominance mechanisms consistent with large mixing angles were 
considered in Refs. [12, [l|]. The dominance of a pseudo-Dirac pair in the solar sector was previously 
displayed in Refs. pq, 12]. Pseudo-Dirac dominance in the atmospheric sector has been recently 
discussed in Refs. 
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